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Abstract Takayasu et al. have revisited the question of stochastic processes with 
muhiphcative noise, which have been studied in several different contexts over the past 
decades. We focus on the regime, found for a generic set of control parameters, in which 
stochastic processes with multiplicative noise produce intermittency of a special kind, char- 
acterized by a power law probability density distribution. We briefly explain the physical 
mechanism leading to a power law pdf and provide a list of references for these results dating 
back from a quarter of century. We explain how the formulation in terms of the character- 
istic function developed by Takayasu et al. [|1| can be extended to exponents /i > 2, which 
explains the "reason of the lucky coincidence". The multidimensional generalization of (|1]) 
and the available results are briefly summarized. The discovery of stretched exponential 
tails in the presence of the cut-off introduced in is explained theoretically. We end by 
briefly listing applications. 
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I. STOCHASTIC MULTIPLICATIVE PROCESSES REPELLED FROM THE 

ORIGIN 



Takayasu et al. recently studied the discrete stochastic equation 

x{t + l) = b{t)x{t) + f{t) , (1) 

as a generic model for generating power law pdf (probability density function). Eq.(|l|) defines 
a stationary process if {Inbit)) < 0. 

In order to get a power law pdf, b{t) must sometimes take values larger than one, cor- 
responding to intermittent amplifications. This is not enough: the presence of the additive 
term f{t) (which can be constant or stochastic) is needed to ensure a "reinjection" to finite 
values, susceptible to the intermittent amplifications. It was thus shown that (0) is only 
one among many convergent ((ln6(t)) < 0) multiplicative processes with repulsion from the 
origin (due to the f{t) term in (|l])) of the form 

x{t + 1) = e^(-W'W*)'/W'-» b{t) xit) , (2) 

such that F for large x{t) (leading to a pure multiplicative process for large x{t)) and 
F (X) for x{t) (repulsion from the origin). F must obey some additional constraint 
such a monotonicity which ensures that no measure is concentrated over a finite interval. 
All these processes share the same power law pdf 

P{x) = Cx-^~^ (3) 

for large x with fi solution of 

{b{tr) = 1 . (4) 

The fundamental reason for the existence of the powerlaw pdf (0) is that In x{t) undergoes 
a random walk with drift to the left and which is repelled from — oo. A simple Boltzmann 
argument shows that the stationary concentration profile is exponential, leading to the 
power law pdf in the x{t) variable. 

2 



These results were proved for the process (|l|) by Kesten using renewal theory and was 
then revisited by several authors in the differing contexts of ARCH processes in econometry 
Q and ID random- field Ising models ^ using Mellin transforms, and more recently using 
extremal properties of the G—Ziarmonzc functions on non-compact groups and the Wiener- 
Hopf technique 0]. Many other results are available, for instance concerning the extremes 
of the process x{t) [0 which shows that x{t) have similar extremal properties as a sequence 
of iid random variables with the same pdf. The subset of times 1 < {te} < t at which x{te) 
exceeds a given threshold xt'' converges in distribution to a compound Poisson process with 
intensity and cluster probabilities that can be explicited 

II. CHARACTERISTIC FUNCTION FOR fi > 2 

Within Renewal theory or Wiener-Hopf technique, the case /x > 2 does not play a special 
role and the previous results apply. In the context of the characteristic function used in [|I| , 
the case > 2 can also be tackled by remarking that the expression of the Laplace transform 
P{f3) of a power law pdf P with exponent /i presents a regular Taylor expansion in powers 
of (3 up to the order / (where / the integer part of fi) followed by a term of the form P^^. Let 
us give some details of this derivation. The Laplace transform 

P{(3) = / dwP{w)e-^'', (5) 
Jo 

applied to @ yields 

= C / dw—— = / dx——. (6) 

We have assumed, without loss of generality, that the power law holds for x > \. Denote / 
the integer part of (/ < /i < / + 1). Integrating by part / times, we get (for C = //) 

+ 7 .) , ^ / dxe~^x'-^. (7) 

(^-l)(^-2)...(^-0^ ^ ^ 
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This last integral is equal to 

rfxe-V-^ = r(/ + l-^)[/?^ + /5'+^7*(/ + l-/i,/5)], (8) 

Jl3 

where F is the Gamma function (r{n + 1) = n\) and 

+OD an 

is the incomplete Gamma function 0. We see that P{P) presents a regular Taylor expansion 
in powers of (3 up to the order /, followed by a term of the form (3^. We can thus write 

P(/?) = 1 + ri/3 + + + r^(3^ + 0(/9'+i), (10) 

with ri = — (x), r2 = ... are the moments of the powerlaw pdf and, reintroducing C, 



where is proportional to the scale parameter C. For small /3, we exponentiate (10) and 
rewrite P{(3) under the form 



P{I3) = exp 
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(11) 



lk=l 

where the coefficient dk can be simply expressed in terms of the r/c's. We recognize in this 
the transformation from the moments to the cumulants. The expression ([TT|) generalizes the 
canonical form of the characteristic function of the stable Levy laws, for arbitrary values of 
/i, and not solely for < 2 for which they are defined. The canonical form is recovered for 
/i < 2 for which the coefficient ^2 is not defined (the variance does not exist) and the only 
analytical term is (for // > 1). This rationalizes "the lucky coincidence" noticed by 

Takayasu et al. [|I| that the results obtained from the characteristic function were found to 
apply numerically for exponents > 2. 



III. MECHANISM FOR THE STRETCHED EXPONENTIAL FOUND BY 

TAKAYASU ET AL. 

To mimick system size limitation, Takayasu et al. introduce a threshold Xc such that for 
x{t)\ > Xc, b{t) < 1 and find a stretched exponential truncating the power law pdf beyond x^. 
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Frisch and Sornette have recently developed a theory of extreme deviations generalizing 
the central limit theorem which, when applied to multiplication of random variables, predicts 
the generic presence of stretched exponential pdf's. Let us briefly summarize the key ideas 
and how it applies to the present context. First, we neglect f{t) in (|1|) for large x{t) {xc is 
supposed much larger than the characteristic scale of f{t)). The problem thus boils down 
to determine the tail of the pdf for a product of random variables. 
Consider the product 

Xn = mim2....m„ (12) 

If we denote p{m) the pdf of the iid random variables mj, then the pdf of X„ is 

Pn{X) ~ [j9(X")]", for X ^oo and n finite. (13) 

Equation ([T3|) has a very intuitive interpretation: the tail of Pn{X) is controlled by the 
realizations where all terms in the product are of the same order; therefore Pn{X) is, to 
leading order, just the product of the n pdf's, each of their arguments being equal to the 
common value X". When p{x) is an exponential, a Gaussian or, more generally, of the 
form oc exp(— Cx''') with 7 > 0, then ( [TB| ) leads to stretched exponentials for large n. For 
example, when p(x) oc exp(— Cx^), then Pn{X) has a tail oc exp(— CnX^/"). 

Expression (|13D is obtained directly by recurrence. Starting from Xn+i = we 
write the equation for the pdf of Xn+i in terms of the pdf's of Xn+i and X„ : 

POO POO 

Pn+li^n+l) = / dXnPn{Xn) / dXn+lp{Xn+l)S{Xn+l — XnXn+l) = 

Jo Jo 



00 







'^^"Pn(X„)pf^) . (14) 



n 

n + l , . , ,^-!- X 



The maximum of the integrand occurs for X„ = {Xn+i)~^ at which Xn = ■ Assuming 
that Pn{Xn) is of the form ([T3|) , the formal application of Laplace's method to (0) then 
directly g 1) is of the same form. Thus, the property (|13D holds for all n 

to leading order in X. See [|ll| for a more detailled derivation. 
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IV. CONCLUDING REMARKS 



The process (|T]) corresponds to a zero-dimensional process. An interesting extension 
consists in taking x to be a function of space (rf-dimension) and time. Qualitatively, we thus 
get a d— continuous infinity of variables x, each of which follows a multiplicative stochastic 
dynamics having the forms (|I|) coupled to nearest neighbors through a diffusion term. Munoz 
and Hwa [jTU[ find numerically a power law decay for the pdf of x in the d-dimensional case. 

Autocatalytic equations lead to multiplicative stochastic equations that are exactly 
tractable in the case of Gaussian multiplicative noise. The process (|T]) also describes 
accumulation and discount in finance, perpetuities in insurance, ARCH processes in econom- 
etry, time evolution of animal population with restocking [§] . The random map (Q) can also 
be applied to problems of population dynamics, epidemics, investment portfolio growth, and 
immigration across national borders . 
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